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Abstract 

The Laser Interferometer Gravftational Observatory (LIGO) has recently reached the end of its 
fifth science run (S5), having collected more than a year worth of data. Analysis of the data is still 
ongoing but a positive detection of gravitational waves, while possible, is not realistically expected 
for most likely sources. This is particularly true for what concerns gravitational waves from known 
pulsars. In fact, even under the most optimistic (and not very realistic) assumption that all the 
pulsar's observed spin-down is due to gravitational waves, the gravitational wave strain at earth from 
all the known isolated pulsars (with the only notable exception of the Crab pulsar) would not be 
strong enough to be detectable by existing detectors. By August 2006, LIGO had produced enough 
data for a coherent integration capable to extract signal from noise that was weaker than the one 
expected from the Crab pulsar's spin-down limit. No signal was detected, but beating the spin-down 
limit is a considerable achievement for the LIGO Scientific Collaboration (LSC). It is customary 
to translate the upper limit on strain from a pulsar into a more astrophysically significant upper 
limit on ellipticity. Once the spin-down limit has been beaten, it is possible to release the constraint 
that all the spin-down is due to gravitational wave emission. A more complete model with diverse 
braking mechanisms can be used to set limits on several astrophysical parameters of the pulsar. 
This paper shows possible values of such parameters for the Crab pulsar given the current limit on 
gravitational waves from this neutron star. 

PACS numbers: 95.55.Ym, 97.60.Gb, 97.60.Jd,04.80.Nn 
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I. INTRODUCTION 

The LIGO observatory has conducted five science runs 
Ul I 01 01 ■ ^ coherent integration to extract signal from 
noise from possible continuous sources, such as lumpy ro- 
tating neutron stars, (using a technique explained in [j|) 
was performed. 

To aid the search it is useful to implement filters matched 
to the signal derived from known neutrons' stars pa- 
rameters such as the celestial coordinates, the rotation 
frequency and the frequency derivatives of the searched 
for source. Known pulsars are very convenient possible 
sources because such parameters are well known from 
radio and X-rays observations. No signal was detected 
from any known pulsars during LIGO's first four data 
runs. The coherent search was limited mainly to sinu- 
soidal signals at frequencies that were twice the rotation 
frequency of the pulsar, typical of a non-axisymmctric 
rotating neutron star. 

Analysis of the data is still ongoing for S5 but a positive 
detection of gravitational waves, while possible, is not re- 
alistically expected from any known pulsar. In fact, it is 
well known that even under the very optimistic and unre- 
alistic assumption that all the spin-down is attributable 
to the emission of gravitational waves, the gravitational 



strain of such radiation reaching the earth would not 
be strong enough to be detectable by existing detectors, 
even after one year of integration time. It is common to 
illustrate this particular fact in a way similar to Fig.l. 
The maximum possible strain from the known pulsars is 
shown as a function of emission frequency and compared 
with the noise level for the an existing LIGO detector and 
planned upgrade. The sensitivity curve is based on the 
reduced noise level achieved after a full year of coherent 
integration (an improvement that is proportional to the 
square root of the integration time) . Almost all the pul- 
sars have maximum possible strains that are below the 
sensitivity curve and therefore are undetectable by the 
present LIGO. The maximum possible strain is an upper 
limit calculated from simple energy conservation argu- 
ments and the value of the observed spin-down of the 
pulsar. The only pulsars with maximum possible strains 
above the curve are the Crab pulsar (PSR B0531-f21) 
and the pulsar J1022-f 1001 that is a member of a binary 
system (the analysis of the gravitational waves from this 
kind of source is more complicated and it will be not con- 
sidered in this letter). 

The spin-down limit for a known pulsar can be calculated 
as in the following. Under the optimistic assumption that 
all the spin-down is due solely to the emission of gravita- 
tional radiation, the pulsar loses energy with a rate given 
by: 
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where H. = 2Tr fr, fr is the rotation frequency, G is the 
gravitational constant, c is the velocity of light, / is the 
moment of inertia and e is the pulsar ellipticity. The radi- 
ated energy depletes the rotational kinetic energy reser- 
voir by an amount E — Qfll. The radiation emitted will 
be in the form of gravitational waves (at twice the ro- 
tation frequency for a simple non-axisymmetric rotator) 
with strain amplitude h that can be determined from the 
time derivative of the energy. In fact, we have that: 
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where / = 2/^ is the gravitational wave frequency and r 
is the distance of the source. 

For the majority of known pulsars the angular frequency 
n, the angular frequency time derivative 51 and the dis- 
tance r are known from radio or X-ray observations. Us- 
ing these values, it is common in the literature to illus- 
trate the spin-down limits as shown in Figure 1. 
The strain spin-down limit for the Crab Pulsar corre- 
sponds to a value of /i < 1.4 x lO"^"*. The spin-down 
limit can be converted into a limit on the ellipticity of 
the pulsar according to the following equation: 
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For the Crab pulsar the spin-down limit on ellipticity is 
e = 7.5 X 10--*. 

The Crab pulsar's emission frequency is very close to 
the power grid typical cycle. Resonances and possible 
frequency up-conversions make the noise floor near the 
Crab's emission frequency higher than what it would be 
without the additional noise due to the power grid inter- 
ference. But enough integration time allowed the LIGO 
detector to reach and beat the spin-down limit. By Au- 
gust 2006 a preliminary upper limit value for the strain 
h = 5.0 X 10~^^ that corresponds to an upper limit on 
the ellipticity e = 2.5 x 10^'' has been claimed by the 
LSC. This result is based on a search that used data 
up to August 25th 2006 when a major glitch occurred 
in the Crab pulsar's timing (creating a possible abrupt 
change of phase) and interrupted momentarily the anal- 
ysis Q . Recently, the results of more accurate and com- 
plete searches were announced [6]. Three searches, with 
different probability methods and physical assumptions, 
were performed giving a range of upper limits on strain 
from 1.7 X 10~^^ to 2.7 x 10~^^, which translates in a 
range of cUipticities values in the range from 9 x 10^^ 
to 1.4 X 10^^. In this paper we mainly use the single 
template limit on ellipticity e = 1.8 x 10~^. This is a 
somehow conservative value because it doesn't depend 
on assumptions made on the orientation of the neutron 
star's spin axis. 
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FIG. 1; The graph shows the gravitational strain hsc (EC 
stands for Energy Conservation or spin-down limit) for the 
known pulsars from the complete ATNF catalog (about 1500 
stars). The two star marks indicate the h upper limit for 
the putative supernova compact remnant SN 1987A (being a 
precessing neutron star it emits a the rotation frequency and 
twice the rotation frequency) Q • The Crab pulsar is above 
the noise curve even for the initial LIGO case. Unfortunately, 
it is very close to the 60 Hz line (caused by resonance with 
the power grid typical cycle.) The quantity hsc was derived 
assuming that all observed spin-down is due to energy loss 
caused by emission of gravitational radiation (and no other 
braking mechanisms.) The strain for the pulsars is compared 
with the noise sensitivity curve (in units of hyHz) for the 
current LIGO detector (actual data for the Hanford detector 
taken during the science run S2 conducted in 2004 and ex- 
trapolated to 1 year of integration time). The expected noise 
curve for Advanced LIGO was computed numerically using 
the free-share code BENCH H. 



II. RATIO OF ENERGY LOSS DUE TO 

GRAVITATIONAL VERSUS 

ELECTROMAGNETIC RADIATION 



A simple model with mixed breaking 
mechanisms 



Once the spin-down limit for a particular pulsar has 
been reached, we can conclude, that in the absence of 
a positive detection of gravitational waves, the star is 
spinning down through a combination of braking mech- 
anisms. For example [9|, we can adopt a model that im- 
plies gravitational and electromagnetic braking. In this 
case the energy loss will be given by: 
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The contribution to the energy loss due to gravitational 
waves is given by: 
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while the contribution to the energy loss due to dipole 
radiation and other electromagnetic processes is given 
by: 
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In general the spin-down can be expressed as a power law 
in terms of rotation frequency fJ as: 



VI = KVf" 



(7) 



where K is the torque function (the quantity 10. is the 
torque acting on the star) and n is the braking index, that 
is telling of the mechanisms at work in spinning down the 
pulsar. For example, in the case of a pure electromagnetic 
dipole n = 3 and n — b for pure gravitational radiation 
(this last result can be easily derived using Eq. (5)). In 
the case of a mixed contribution of braking mechanisms, 
n will be different from the above values. In this case we 
will have: 
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where Kf^m and riem are the electromagnetic torque func- 
tion and braking index and Kqw and new are the equiv- 
alent gravitational quantities. In principle, the power law 
n can be calculated from observational parameters of the 
pulsar. In fact, using Eq. (4), we can derive the expres- 
sion: 
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Just for few pulsars (about four) we can measure the 
value of VI reliably, therefore a calculation of n from ob- 
servable quantities is usually unpractical. 
The Crab pulsar is one of such pulsars for which the value 
of n is measured and it is equal to 2.51. This allows us 
to set some interesting upper limits on the energetics of 
the pulsar and on the index of the braking mechanisms 
at work. 

Following Q, we can define the ratio Y{Vl) between the 
energy loss due to gravitational waves and the energy loss 
due to electromagnetic processes as: 
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The braking index n can be expressed in terms of the 
quantity Y as: 



l + Y ' 
and then solving for Y as a function of n we have: 



Y{n) = 
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Using Eq. (5), (10) and (12) we can calculate an upper 
limit on Y{n) from the upper limit on ellipticity using 
the following expression: 
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0.058 


upper limit 


^eni 


2.36 


lower limit 


Eqw (in percentage) 


5.55 


upper limit 


-Bern (in percentage) 


94.44 


lower limit 



TABLE I: Upper and lower limits of energy ratio U„ be- 
tween gravitational radiation and electromagnetic braking, 
electromagnetic braking index n^rn and upper and lower limits 
of gravitational and electromagnetic radiation contributions 
Egw and Eem to the total energy loss. 



then solving for Y{n): 



Y{n) 



I- A 
where the expression A is equal to: 
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Finally, we can calculate a lower limit on the electromag- 
netic braking index Uf^rn using Eq. (12): 
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We can also write an expression for the ratio of the 
energy contribution of the electromagnetic and gravita- 
tional braking to the total energy loss Etot'- 
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Using the Crab pulsar observational parameters n = 
2.51, VI = 186.96 s-i, VI = -2.33 x lO^^s"^ ^^d latest 
ellipticity limit e — 1.8 x 10^^ from the LSC, we obtain 
current values for the upper and lower limits for the ratio 
Ym the electromagnetic braking index riem and an upper 
limit on the contribution of gravitational radiation to the 
total energy loss as illustrated in Table 1. 



B. Energetics of the Crab pulsar's spin-down 

A large body of literature discusses the topic of the en- 
ergy budget of the Crab pulsar and the surrounding Crab 
Nebula. See Ref. [i3|, [iil] for comprehensive reviews. 
The wind nebula acts as an effective calorimeter that al- 
lows us to measure the energy interactions between the 
pulsar and its environment. The mechanical spin-down is 
considered the main power source for both the pulsar and 
the Nebula luminosity. This power is about 4.77 x 10^^ 
ergs s~i. The pulsar total luminosity in different elec- 
tromagnetic wavelengths accounts just for a few percent 



of the total observed spin-down. For example, the pulsar 
radio luminosity is about 6.9 x 10'^^ ergs [I^l or 1.5x10^^ 
of the mechanical loss of energy, while the X-ray lumi- 
nosity is about 10~^. The total luminosity of the Nebula 
is about 1.3 x 10^^ ergs s~^. It is well established that 
the main mechanism that allows the pulsar to transfer 
a large part of its mechanical energy to the Nebula is 
synchrotron radiation. For this process an efhciency of 
about 30 percent is sufficient to explain the observed to- 
tal Nebula's luminosity. The remaining 70 percent of 
the energy loss is somehow not showing up in either the 
Nebula or the pulsar electromagnetic spectrum. In the 
mixed braking model that was discussed in the previous 
section, the upper limit on the efficiency of the gravita- 
tional braking was about 6 percent. Therefore, this value 
of the efficiency is consistent with the Crab pulsar's and 
Nebula's energy budget as based on current observations. 
Instead, the assumption usually adopted to calculate the 
spin-down limit (all the energy loss is due only to grav- 
itational waves) is, of course, in flagrant contrast with 
the above arguments. In fact, a large amount of energy 
is released into the Nebula. This energy cannot be in the 
form of gravitational waves that do not interact readily 
with matter. 

The Crab pulsar's energy loss due to gravitational waves, 
implied by the most recent upper limit from LIGO, is 
now more consistent with other astrophysical observa- 
tions, as the pulsar and Nebula electromagnetic spec- 
trum. The implied upper limit on ellipticity itself is still 
too large from a theoretical point of view. The current 
value e = 1.8 X 10^* is at best at the limit of what the 
most exotic neutron star models allow [l3|. 
As better limits on the gravitational strain of the Crab 
pulsar will be set by future LIGO data, we will have 
tighter limits on the efficiency of the gravitational wave 
contribution to the spin-down. At the same time, the 
contribution of the gravitational braking to the overall 
energy loss will become quickly less quantitatively sig- 
nificant as the efficiency depends on the square of the 
ellipticity. An improvement of a factor of 10 in strain 
sensitivity (easily reachable by Advanced LIGO) would 
mean an improvement on the upper limit on ellipticity of 
the same order. In this case, the upper limit on the effi- 
ciency of the gravitational wave braking will be about 100 
times smaller than the current one. At that point, the 
gravitational wave luminosity upper limit will be compa- 
rable with the pulsar's X-ray luminosity (and yet about 
ten thousand times higher than the radio luminosity). 
See Fig. 2. 



III. A MORE COMPLETE MODEL WITH n < 3 

The dipole radiation model is the most traditional 
braking mechanism suggested [1J| and one of the most 
successful in explaining many of the pulsars' observed 
properties. 
An interesting alternative model is proposed by Chan 
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FIG. 2: Efficiency of the Crab pulsar's gravitational radiation 
braking versus ellipticity. The ellipticity is in the range be- 
tween the upper limit given by the current LIGO observations 
and possible limits given by future observations. The differ- 
ent curves represent different values of the moment of inertia. 
7o = 1 X 10 '^ gem is the canonical value for the moment of 
inertia. The ellipticity and a neutron star's moment of inertia 
are dependent parameters given a certain observational upper 
limit on the gravitational strain. See section III. 



and Li (2006) ^M,- The authors assume in their model a 
conventional dipole torque and, to make the overall brak- 
ing index consistent with a value less than 3, a fallback 
disk torque is invoked. It is possible to add a gravita- 
tional torque produced by an asymmetry in the crust of 
the Crab pulsar as given by the LIGO ellipticity upper 
limit. In this case, the energy loss can be described by 
the following equation: 
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where the energy loss contribution due to the fall back 
disk can be obtained from: 
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where c is the speed of light, R is the radius of the neutron 
star, t is the historical age of the star, f3 is the pressure 
scale height to radius ratio, 7 is the light-cylinder radius 
of the pulsar and a is the viscosity parameter. As in [l5| 
we use, for our calculations in this letter, the following 
values for these parameters: (3 = 0.35,0.45,0.55, 7 = 1 
and a — 0.001. The parameter M is the mass inflow 
rate (assumed to be constant). Once the other param- 
eters are set we are going to calculate lower limits on 
the parameter M for a given value of the free parameter 
/3. Using equation (15) it is possible to show that the 
braking index n can be rewritten as in the following: 
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FIG. 3: The fallback disk mass inflow rate M as a function 
of th moment of inertia / of the pulsar and the dimensionless 
parameter /3. 



where r = — f2/(2fi) is the canonical age of the pulsar, 
^GW — ^gvk/(^^) is the efficiency of the gravitational 
torque and r\T^ — (lTd/{I^) is the efficiency of the disk 
braking torque. Equation (15) and (17) can be solved as 
a system to find unique values for the lower limits on the 
magnetic field B of the pulsar and the fallback disk mass 
inflow rate M given the parameter /? and the upper limit 
on ellipticity of the Crab pulsar derived from the LIGO 
null observation. 

To further our analysis it is important to point out that 
in actuality the ellipticity e and the moment of inertia / 
are not independent parameters in our model. In fact, 
the gravitational strain upper limit sets an equivalent 
upper limit only on the quadrupole moment le. It is 
then customary to use a canonical value of / = 10**^ g 
cm^ to calculate an upper limit on ellipticity. But the 
possible values for /, suggested by theory or observation, 
span a range within a factor of few from the canonical 
value [lq |. Therefore, we express our results in terms of 
dependency on the moment of inertia /. Figure 3 shows 
the mass inflow M as a function of the moment of inertia 
/ for different realistic values of the parameter /3. For 
example, for a value of /3 = 0.035 and a canonical value 
of / we have B = 1.07 x 10^^ gauss and M = 2.87 x 10^^ 
gs~l. 

Figure 4 compares the efficiencies of the different brak- 
ing mechanisms as a function of the angular momentum 
/. The greatest efficiencies for the gravitational wave 
braking is given by small values of /. Notice that there is 
a limiting value of about / = 2.2 x 10'''* g cm^ that still 
allows for the presence of a gravitational wave torque and 
the existence of a fallback disk. Values of / smaller than 
this give an unrealistic negative mass inflow. 
These are just examples of how upper limits based on 
gravitational wave observations can start to constrain as- 
trophysical models of neutron stars. For what concerns 
the energetics of this model, using the canonical value of 



FIG. 4: The efficiency of the different breaking mechanisms 
for the fallback disk model as a function of the pulsar moment 
of inertia I. Small moments of inertia correspond to large 
efficiency of the gravitational torque. The limiting value of 
I = 2.2 X 10*'* g cm^ corresponds to the smallest value of / 
consistent with the presence of a fallback torque. 



the moment of inertia, we have that the electromagnetic 
braking efficiency is equal to about 80 percent, the grav- 
itational braking efficiency is 6 percent and the efficiency 
of the disk braking torque is about 14 percent. Given 
the large ellipticity assumed in this particular model the 
gravitational torque is actually comparable in magnitude 
to the fallback disk torque. 



IV. CONCLUSIONS 

The reaching and beating of the Crab pulsar's spin- 
down upper limit is an important milestone for the LIGO 
project. The current strain sensitivity upper limit can 
be converted to give an upper limit for the ellipticity of 
1.8 X 10-4. 

This upper limit for ellipticity is still too large to be re- 
alistic from an astrophysical point of view. However, we 
can start to make interesting astrophysical statements on 
the braking mechanisms and energetics at play in slowing 
down the star. 

In fact, we can now free the constraint that all the en- 
ergy loss is due to the emission of gravitational radiation 
and assume a mixed braking mechanism with contribu- 
tions from different physical phenomena (gravitational 
radiation, electromagnetic dipole radiation, winds, mag- 
netospheric currents, fallback disk torques, etc). 
Different models can be adopted. One of the models we 
have studied assumes that the electromagnetic braking 
index is different from the conventional n = 3. We can 
use the upper limit on the ellipticity to calculate the con- 
tribution to the energy loss due to gravitational radiation 
and calculate the parameter Uem, the braking index due 
to electromagnetic processes. Our current result for the 



Crab pulsar is that the ratio between the energy loss due 
to gravitational radiation to the energy loss due to elec- 
tromagnetism is 0.058. This implies that the contribution 
of the gravitational radiation to the total energy loss has 
to be less than 6 percent. Finally, the electromagnetic 
braking index is 2.36 (this is a lower limit). 
Alternatively, we can use a model that allows for a clas- 
sical electromagnetic dipole but assumes other braking 
mechanisms to account for a overall breaking index dif- 
ferent from n = 3. In this case we can show that the 



efficiency of the gravitational wave braking (given the 
current limits on ellipticity from LIGO) is comparable to 
that of a fallback disk torque with realistic parameters. 
The limit on ellipticity can be used to set limits on the 
mass inflow for the fall back disk. It would be neces- 
sary to have a mass inflow M of about 2.87 x 10^* 5s^l 
to make the model consistent with the limit on elliptic- 
ity. Even more interesting conclusions on the limits of 
the braking mechanisms and energetics of pulsars will be 
reached with advanced detectors. 
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